Phase transitions detached from stationary points of the energy landscape 
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The stationary points of the potential energy function V are studied for the </> 4 model on a 
two-dimensional square lattice with nearest-neighbor interactions. On the basis of analytical and 
numerical results, we explore the relation of stationary points to the occurrence of thermodynamic 
phase transitions. We find that the phase transition potential energy of the 4 model does in 
general not coincide with the potential energy of any of the stationary points of V. This disproves 
earlier, allegedly rigorous, claims in the literature on necessary conditions for the existence of phase 
transitions. Moreover, we find evidence that the indices of stationary points scale extensively with 
the system size, and therefore the index density can be used to characterize features of the energy 
landscape in the infinite-system limit. We conclude that the finite-system stationary points provide 
one possible mechanism of how a phase transition can arise, but not the only one. 
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The stationary points of the potential energy function 
or other classical energy functions can be employed to 
calculate or estimate physical quantities. Well-known ex- 
amples include transition state theory or Kramers's reac- 
tion rate theory for the thermally activated escape from 
metastable states, where the barrier height (correspond- 
ing to the difference between potential energies at certain 
stationary points of the potential energy function) plays 
an essential role. More recently, a large variety of related 
techniques has become popular under the name of energy 
landscape methods with applications to many-body 
systems as diverse as metallic clusters, or biomolecules 
and their folding transitions. While the mentioned ap- 
plications focus mostly on the numerical investigation of 
finite systems, the analysis of stationary points has also 
proved useful for analytical studies of TV-body systems 
in the thermodynamic limit. One field of research where 
such methods have been fruitfully applied is disordered 
systems undergoing a dynamical glass transition 

Another line of research based on stationary points but 
focusing on equilibrium phase transitions in the thermo- 
dynamic limit N — > oo, dates back to about the same 
time [3j]. This approach, originally formulated in terms 
of topology changes of configuration space submanifolds, 
can be rephrased in terms of stationary points of the po- 
tential energy function V, i.e. configuration space points 
q s satisfying W(q s ) — 0. The underlying idea can be 
understood as follows Thermodynamic equilibrium 
properties are encoded in the thermodynamic limit value 
of the microcanonical configurational entropy 

s N (v) = lln fc[V( X )-Nv]d X = iln £ (1) 

where V denotes configuration space and dx its vol- 
ume measure, T, v c T is the hypersurface of constant- 



potential energy V = Nv, and dS stands for the (N — 1)- 
dimensional Hausdorff measure on £„. At a stationary 
point, we have W = 0, the integrand on the righthand 
side of Q]) diverges, and we may expect the stationary 
point to give an important contribution to the integral. 
Indeed, it has been shown that, for finite N, every sta- 
tionary point q s of V induces nonanalytic behavior in 
sn{v) precisely at the potential energy of the stationary 
point, v = V(q s )/N [|. 

Nonanalyticities of thermodynamic functions are hall- 
marks of phase transitions. Having observed that sta- 
tionary points of V cause nonanalyticities in the finite- 
system entropy sjv 5 it seems natural to ask whether they 
may also be responsible for nonanalytic behavior in the 
infinite-system entropy, i.e., for the occurrence of phase 
transitions. While for finite N every stationary point 
of V induces a nonanalyticity in sjv, the majority of 
the nonanalyticities does not survive the thermodynamic 
limit. Two questions turn out to be central to the un- 
derstanding of these observations. (1) Under what con- 
ditions can the nonanalytic behavior of sn , induced by a 
stationary point of V, survive the thermodynamic limit? 
A possible scenario has been depicted in [H, El, and it 
turns out that the Hessian determinant of the potential 
energy function, evaluated at the stationary points, is 
crucial for discriminating whether or not the stationary 
points can induce a phase transition in the thermody- 
namic limit. For some models this insight has proved 
particularly useful in that it facilitates the analytic com- 
putation of phase transition energies even in the absence 
of an exact thermodynamic solution Q. (2) Are station- 
ary points necessary for a phase transition to take place? 
Or is there another mechanism, distinct from the one 
sketched above, which can give rise to a phase transition? 
This question was addressed by a theorem in Q, claim- 
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ing that, for a large class of systems with short-range 
interactions, the presence of stationary point is indeed 
necessary for a phase transition to occur. 

The model study presented in this Letter puts this al- 
leged relation of stationary points of the potential energy 
function V and thermodynamic phase transitions to the 
test. We apply three different methods to extract infor- 
mation about the stationary points of the potential en- 
ergy function of the two-dimensional nearest-neighbor 4 
model. This model is known to have a continuous phase 
transition which is in the universality class of the two- 
dimensional Ising model but, in contrast to that model, 
is amenable to energy landscape methods by virtue of 
its continuous configuration space. The results from two 
complementary numerical techniques, as well as a rig- 
orous analytical upper bound, all provide evidence that 
stationary points occur only at nonpositive potential en- 
ergies. 

The implications of this finding are significant. Among 
other things, they show by counterexample that the the- 
orem on the relation between stationary points and ther- 
modynamic phase transitions, allegedly proven in Q, 
is incorrect. Instead, we observe that, even for short- 
range interacting systems, thermodynamic phase transi- 
tions can occur at energies not related to any station- 
ary points. As a consequence, further mechanisms of 
how phase transitions arise, possibly not related to finite- 
system stationary points, must exist, and we will com- 
ment on possible scenarios towards the end of this Letter. 

Two-dimensional nearest-neighbor <fi model. — On a fi- 
nite square lattice A C 1? consisting of N = L 2 sites, a 
real degree of freedom fa is assigned to each lattice site 
i e A. By Af(i) we denote the subset of A consisting 
of the four nearest-neighboring sites of i on the lattice 
under the assumption of periodic boundary conditions. 
The potential energy function of this model is given by 
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where q — (qi,.. • ,.<?at) denotes a point in configuration 
space L = R N [9]. J > determines the coupling 
strength between nearest-neighboring sites and the pa- 
rameters A, fi > characterize a local double-well poten- 
tial experienced by each degree of freedom. 

In the thermodynamic limit N — > oo this model is 
known to undergo, at some critical temperature T c , a 
continuous phase transition, in the sense that the con- 
figurational canonical free energy f(T) is nonanalytic at 
T = T c . The transition is from a "ferromagnetic" phase 
with nonzero average particle displacement to a "param- 
agnetic" phase with vanishing average displacement (see 
10] for more details as well as for Monte Carlo results). 

Instead of T c , it is more adequate for our purposes to 
compare to the critical potential energy per lattice site 
v c of the transition Both quantities are unambigu- 
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FIG. 1. (Color online) Monte Carlo results for the two- 
dimensional (f> 4 model Q with A = 3/5, fj 2 = 2, and J — I. 
System sizes N = L x L are plotted with L ranging from 32 to 
128. Left panel: the canonical average (v) of the potential en- 
ergy per lattice site as a function of the temperature T. Right 
panel: The canonical specific heat c = ((V 2 ) - (V) 2 )/(NT 2 ), 
plotted as a function of (v) . The lines between the data points 
serve as a guide to the eye. 



ously related to each other in the thermodynamic limit 
via the caloric curve v(T), independently of the statisti- 
cal ensemble used. The critical potential energy v c is less 
frequently studied, in fact the only data we could find in 
the literature are from Monte Carlo simulations of fairly 
small system sizes N = 20 x 20 in [l§], with parameter 
values A = 3/5, p? = 2, and J = 1. We mostly use the 
same parameter values in the following. Since the value 
of v c is a crucial benchmark when relating our stationary 
point analysis to the phase transition of the cf) 4 model, we 
have performed standard Metropolis Monte Carlo simu- 
lations for somewhat larger system sizes up to 128 x 128 
and 10 7 lattice sweeps. We have sampled, among other 
observables, the canonical average (v) of the potential en- 
ergy per particle, and plots of some of the Monte Carlo 
results are shown in Fig. [TJ The quality of the data is 
sufficient to establish, beyond any reasonable doubt, a 
critical potential energy v c 2.2 well above zero. 

Numerical continuation method. — For J = 0, i.e., in 
the absence of coupling, the stationary points q s of the 
potential 1^ in (|2]) can be calculated analytically without 
difficulty, obtaining the 3 N solutions q s — (q\, . . . ,q s N ) 
with q S j £ {0, ±-\/6/i 2 /A}. Knowledge of these solutions 
permits us to continue them to J > by numerical con- 
tinuation; see [13} for a description of the homotopy con- 
tinuation method we have actually been using. Under 
certain conditions on the initial (decoupled) and final 
(coupled) potentials, this method is known to yield all 
stationary points of V. However, since the number of 
stationary points for J = grows exponentially with the 
number of lattice sites, such an analysis is restricted to 
fairly small system sizes. 

We have used the homotopy continuation method to 
compute all stationary points of V for various values of 
J and system sizes up to 4 x 4, and the results show 
the following features. First, upon increasing the cou- 
pling constant J, the number #(<? s ) of stationary points 
dramatically decreases from 3^ for J = to only 3 sta- 
tionary points for larger J. This behavior is illustrated 
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FIG. 2. (Color online) Numerical results from the homotopy 
continuation method. Left panel: The number of stationary 
points of V for N = 3x3, plotted logarithmically as a function 
of the coupling J. Right panel: For all 10 288 973 stationary 
points q s of a 4 x 4 lattice with J = 0.2, the scaled Hessian 
determinant D = | defHy(q s )| 1//JV is shown vs the stationary 
value v B = V(q B )/N, demonstrating that v B < for all q s . 



for N — 3 x 3 in the left panel of Fig. [3] Second, the 
stationary values v s = V(q s )/N, i.e., the potential en- 
ergy per lattice site evaluated at a stationary point, is 
found to be never larger than zero. This is illustrated for 
N = 4 x 4 in the right panel of Fig. [2 

Newton-Raphson method. — We use Monte Carlo dy- 
namics in configuration space to generate a large set of 
initial states, and apply the Newton-Raphson method to 
these initial states to find stationary point of V. In con- 
trast to the homotopy continuation method, there is in 
general no way of knowing whether all stationary points 
have been found. An advantage, however, is that the 
Newton-Raphson method can be applied to system sizes 
much larger than 4x4. 

We have used the routine newt from a globally 
convergent version of the Newton-Raphson method, to 
compute stationary points of (/> 4 lattices of sizes up to 
32 x 32. For small J, the number of stationary points 
for such system sizes is huge and only a tiny fraction of 
them can be tracked down. For sufficiently large J [IH, 
however, only three stationary points are found: the two 
global minima q s = (q\, . . . ,q s N ) where all q* = ^J&JjF/X, 
respectively — -\/6/j 2 /A, and a stationary point of index 
1 at qj — for all j. Moreover, we again find that all 
stationary values obey v s < 0. 

Analytical upper bound on the stationary values. — The 
stationary points of the potential V are the solutions of 

^L = Ll + (4J-» 2 )q k -J E *=0- (3) 

Although it is not feasible to solve this set of N cou- 
pled nonlinear equations explicitly, the potential energy 
at a stationary point can be determined by rewriting the 
potential ^ in the form 
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energy at a stationary point q s = (q±, . . . , gjv), 
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Since A > 0, the potential energy per lattice site at any 
stationary point is bounded above by zero, 



V s = V(q s )/N < 0. 



(6) 



Comparison with earlier results. — Our findings, and, 
in particular, the fact that the stationary values v s are 
nonpositive, disprove earlier results on the relation be- 
tween thermodynamic phase transitions and stationary 
points of V. These earlier results were phrased in terms 
of topology changes of certain submanifolds in configu- 
ration space, but with the help of Morse theory we can 
rephrase all statements in terms of stationary points. 

In 2004, Franzosi and Pcttini announced, and allegedly 
proved under some conditions on the potential V, a nec- 
essary condition for a thermodynamic phase transition 
to occur [§]. In the language of stationary points, their 
claim can be phrased as follows: If there exists an interval 
[a, b] such that, for all system sizes N larger than some 
constant Nq, the stationary values v s = V(q s )/N corre- 
sponding to the stationary points q s of V all lie outside 
that interval, then in the thermodynamic limit neither a 
first- nor second-order thermodynamic phase transition 
can occur at critical potential energies v c G (a, b). 

In short, stationary points in the vicinity of some v c 
are claimed to be necessary for a phase transition at v c . 
The nearest-neighbor 4 potential ([2]) satisfies all require- 
ments on V demanded by this theorem. Accordingly, 
based on the fact that all stationary values are nonpos- 
itive (H]), the theorem asserts that the critical potential 
energy of the second-order phase transition of the model 
cannot be positive. This prediction is in contradiction 
to the thermodynamic properties of the model and the 
theorem in 0] is falsified by means of a counterexample. 

A numerical study of the configuration space topol- 
ogy of the two-dimensional nearest- neighbor </> 4 model 



was published by Franzosi et al. in 16|. The authors 
reported results for the Euler characteristic \ ( a topo- 
logical invariant) of the constant-potential energy shells 
in configuration space, finding a pronounced kink of \ 
as a function of v in the vicinity of the transition energy 
(Fig. 3 of [16]). From the absence of stationary points 
at positive energies, Morse theory allows one to conclude 
that x( v ) 1S rigorously constant for v > 0. The kink 
observed in Fig. 3 of [16( is therefore an artifact of the 
numerical method used. 

More on the shape of the energy landscape. — Contrary 



to the claims in [8|, |16J, we have seen that the energy 



Then, substituting © into (jj), we obtain the potential 



landscape in the vicinity of the critical potential energy 
v c of a phase transition can be locally trivial, i.e., free of 
stationary points with potential energies v s in the vicin- 
ity of v c . For v < and large system sizes, however, 
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FIG. 3. (Color online) Left panel: The relative index i of the 
stationary point q s = (0, . . . , 0), plotted as a function of the 
linear system size L = y~N for various couplings J. From top 
to bottom: J = 0.1, 0.36, 0.5, 0.64, 1, 2, 10. Right panel: 
The large-system limit ioo of the relative index, plotted as a 
function of the coupling constant J. 

when stationary points abound, it is more difficult to ex- 
plore the shape of the energy landscape. As a first step 
towards this aim, we study the properties of the station- 
ary point q s — (0,...,0), corresponding to the largest 
stationary value v s = 0. Its index /, i.e., the number of 
negative eigenvalues of the Hessian matrix of V evaluated 
at q s , characterizes the change of the constant-energy 
shell around v — 0. With increasing system size N, we 
observe that the relative index i = I/N converges to a 
finite value, and this value i^ = limjv_ ! . 00 i/N depends 
on the coupling J (see Fig. [3]). The existence of such 
a limiting value is a good starting point for discussing 
the properties of the potential energy landscape in the 
infinite-system limit. In fact, we can deduce that the 
potential energy landscape does not approach a simple 
iV-dimensional generalization of a double-well potential 
in the large-system limit, as one might naively have ex- 
pected. Such a double well, having two degenerate min- 
ima at the ground state energy and a stationary point 
of index I = 1 at v = 0, would yield ioo = for the 
stationary point q s = (0, . . . , 0). 

Conclusions. — By analytical and numerical methods, 
we have probed certain features of the potential en- 
ergy landscape of the two-dimensional (j) 4 lattice model. 
The model's phase transition was found to occur at en- 
ergies well separated from the stationary values v s = 
V(q s )/N of the potential V (or, equivalently, from topol- 
ogy changes in configuration space). These findings fal- 
sify a theorem put forward in [8j which claims that sta- 
tionary points q s with stationary values v s — v c are neces- 
sary for a phase transition to occur at a critical potential 
energy v c . Since our results imply that the constant- 
potential energy shells are simply-connected for v > 0, 
the symmetry breaking phase transition is found to be 
driven by a concentration-of-measure effect, but is not re- 
lated to the connectivity of the underlying finite-system 
energy shells. In other words, even if the integral on the 
right-hand side of ([T]) depends smoothly on v, the limit- 
ing procedure N —¥ oo can, contrary to the claim in Q, 
destroy this smoothness. 

Accordingly, we conclude that the finite-system sta- 



tionary points provide one possible mechanism of how a 
phase transition can arise, but not the only one. If the 
stationary points are at the basis of the transition, the 
earlier mentioned criterion Q based on the Hessian de- 
terminant at stationary points remains valid and can be 
applied to analyze the phase transition. However, other 
scenarios are possible and open up interesting perspec- 
tives. One possibility is that, even if V has no stationary 
points when considered on configuration space, it might 
have such points when treated as a function of N com- 
plex variables. If some of these stationary points, in the 
thermodynamic limit, approach the real axis and satisfy 
the Hessian determinant criterion at the same time, they 
should be capable of inducing a phase transition, despite 
the absence of stationary points in the real (noncomplcx) 
configuration space. This and other alternative scenarios 
could open up possibilities for deriving new criteria on 
the existence or absence of phase transitions, as well as 
for analytic methods for computing transition energies. 
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